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; 1 Introduction 

We examine the impact of time-dependent mass injection and heating on the evolu- 
tion of the interstellar medium (ISM) in elliptical galaxies. As the large and luminous 
ellipticals have supermassive black holes at their cores, which were probably much less 
massive in the young universe, feeding these black holes is essential. The evolution of the 
| ISM, controlled by supernov a and/or AGN heating as well as cooling by bremsstrahlung 

■ (jMathews fc Brighenti 2003h . could provide an abundant matter supply for this if a suf- 



ficient fraction of the ISM flows towards the center. 
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l/-) ( The hot phase of the ISM has typical temperatures around 10 K at the galactic 

cores. While stellar winds and planetary nebulae are widely seen as the sources of the 
■ interstellar gas, there is no generally established model to explain the heating of the ISM 

(supernovae, AGN heating, collisions). We concentrate on supernovae as a source of the 
O ' heating. 

CO 

2 Model and Basic Equations 

t > ■ 

We use simple Plummer models as well as combined Hernquist/NFW models for the 
mass distribution of our galaxies, with p*(r) being the stellar mass density. In addition, 
we have a black hole with mass Mh at the center. Most calculations are performed in 
one spatial dimension, some in 2D including r andom disturbances, w ith the magneto- 
hydrodynamic code NIRVANA by Udo Ziegler ( Ziegler k, Yorke 1997^ . 



^ + V • (pv) = ap* (1) 
+ V- (pv®v) = -Vp-pV<S> (2) 

at 

— + V ■ (ev) = -pV ■ v + ap*c v T - p 2 A (3) 
at 
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Figure 1. Typical temperature profiles of steady state outflow and inflow solutions. The two types 
of solutions can be easily distinguished in this figure: While inflows generally have low temper- 
atures and positive gradients, galactic winds show high temperatures and negative gradients, 
which might be very difficult to observe directly as the densities are much lower, a = 10~ 19 s _1 , 
inflows: T = 5 • f 5 K . . . 10 7 K, outflows: T = 1 • 10 s K . . . 9 ■ 10 s K. 

Here p, u, T, p and e (energy density) are the state variables of the gas, p 2 A is the 
cooling rate, and c v To is the supplied energy relative to the injected matter ap+, so the 
heating is controlled by the value of To. The mass injection by stellar mass loss as well 
as the supernova heating are assumed to be proportional to the stellar density because 
of their physical origins. 

3 ISM Simulations 

For constant mass injection and heating we find steady state solutions, which can be 
characterized as inflow and outflow types (Figure 1). Inflow solutions show temperatures 
near the virial temperature and subsonic gas flows towards the center (cooling flows) 
independent of the heating rate. Outflows on the other hand generally are supersonic, 
of high temperature and low gas density. Their velocity and temperature profile is 
independent of the mass injection rate. 

Although there is a transition region between inflow and outflow in the (a, To) para- 
meter space where we have spatially bimodal flows, this region seems to be narrow. If 
there exists at least a partial inflow, then a significant fraction of the globally injected 
mass flows towards the center. 

Figure 2 shows the distribution of the two solution types in the parameter space. 
There obviously exists a critical heating parameter To iCr it which separates these types. 
This can be explained easily by checking the energy balance. The ISM evolution in 
our simulations with time-dependent mass injection and heating can be understood by 
examining the evolution of the model in the (a, To) parameter space. 

The time-dependence of mass injection follows from the initial mass function of the 
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Figure 2. Distribution of inflow (filled symbols) and outflow (open symbols) solutions in the 
(a, To) parameter space for our model galaxy (dark matter contribution is not included). The 
mass injection parameter a directly corresponds to the overall mass injection rate in the galaxy, 
while the heating parameter To depends on the supernova rate (in SNu) and the mass injection. 
1 SNu corresponds to 1 supernova per 100 yr and 10 10 L Q . The diagonal lines represent constant 
supernova rates. 
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Figure 3. Evolution of SN la heating in the (a, To) parameter space for a powerlaw behavior 
at late times in the galaxy: supernova rate oc t~ s and mass injection a oc t . The arrows 
indicate the a-position of the galaxy at a given time. 
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stellar population generated by a monolithic star burst. The amount of m ass injection 
seems to be in good agreement with observations ( Breeman fc Athev 2004h . 

The time-dependence of the supernova rate is rather unknown. During the initial 
starburst there is much contribution from SNe II, but soon after SNe la are the only ones 
occurring. At late times, in passively evolving galaxies, the evolution of the SN la heating 
rate is often assumed to be a powerlaw (oc t~ s ), but observational confirmation for this 
is still missing so that the powerlaw exponent could vary over a large range. Figure 3 
shows the evolutionary tracks of the model for three different powerlaw exponents. 

The critical heating paramete r is governed by the vir ial mass of the galaxy, and su- 
pernova rates in local ellipticals ( Cappellaro et al. 19991 ) suggest that our model galaxy 
of lO n M virial mass should be near the turnover region between inflow and outflow. 
For early times, we find a galactic wind for s = 1.7 and inflow for s = 0.9. This shows 
that higher supernova rates in the past (s > 0) are not sufficient to drive galactic wind s 
- they have to increase stronger (s > 1.3) than the mass injection ( Ciotti et al. 199lh . 

For a galaxy with much higher or much lower mass, the position of the turnover 
region is much more important than the evolution of the supernova rate so that high- 
mass galaxies would generate inflows with high X-ray luminosities while low-mass dwarf 
galaxies would have strong galactic winds and unobservable X-ray luminosities. 



4 Impact of the Initial Starburst 

At very early times in the galaxy, during its formation, the abovementioned powerlaw 
behavior of course breaks down. Figure 4 shows the evolution of the parameters a and To 
during galaxy formati on based on two very simple scenarios with data from Starburst99 
( Leitherer et al. 1999^ : The galaxy is instantaneously formed from a gas cloud or the 
galaxy is formed during 500 Myr by a constant star formation rate. The mass injection 
a decreases with increasing formation timescale, but the heating parameter To remains 
the same because To is proportional to the heating rate per injected mass. 

The extremely high supernova rates during the starburst lead to high heating rates, 
but this does not necessarily lead to galactic winds. Figure 4 shows that To ~ 10 8 K 
during the starburst. This would drive a galactic wind for the model galaxy, but couldn't 
do so for a galaxy with a virial mass of > 10 12 M as the threshold-T) for winds increases 
proportional to the galaxy mass. The reason is that a huge amount of injected matter 
would have to be lifted out of the galaxy's gravitational potential. 



5 Effects on Black Hole Growth 

The evolution of the ISM might be a crucial factor for fuelling central black holes. If 
a galactic wind arises and stays for a long time, this could effectively avoid black hole 
growth because of low gas densities and high temperatures. A big galaxy in inflow state, 
on the other hand, could easily accrete some 10 10 M because the total injected gas mass 
is of the order of the stellar mass of the galaxy. 
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Figure 4. Evolution of the mass injection and heating parameters a and To during galaxy for- 
mation for different metallicities Z, derived from Starburst99. The time-dependence of mass 
injection (top) and heating parameter (bottom) are shown for instantaneous formation of the 
galaxy and continuous formation during 500 Myr. Heating by SNe la is not included. 



Depending on the galaxy mass and the evolution of the heating rate, there is a rich 
variety of possible scenarios for the interstellar gas. Inflow at early times can produce 
supermassive black holes very quickly and possibly resulting activity could provide a 
feedback through additional heating which could turn the inflow into an outflow. The 
impact of further star formation in cooling flows could be a similar feedback mechanism 
and induce a cyclic behavior. 
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